Charged-phonon theory and Fano effect in the optical spectroscopy of bilayer graphene 
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Since their discovery, graphene-based systems represent an exceptional playground to explore the 
emergence of peculiar quantum effects. The present paper focuses on the anomalous appearence 
of strong infrared phonon resonances in the optical spectroscopy of bilayer graphene and on their 
pronounced Fano-like asymmetry, both tunable in gated devices. By developing a full microscopic 
many-body approach for the optical phonon response we explain how both effects can be quan- 
titatively accounted for by the quantum interference of electronic and phononic excitations. We 
show that the phonon modes borrow a large dipole intensity from the electronic background, the 
so-called charged-phonon effect, and at the same time interfer with it, leading to a typical Fano 
response. Our approach allows one to disentangle the correct selection rules that control the rel- 
ative importance of the two (symmetric and antisymmetric) relevant phonon modes for different 
values of the doping and/or of the gap in bilayer graphene. Finally, we discuss the extension of the 
same theoretical scheme to the Raman spectroscopy, to explain the lack of the same features on the 
Raman phononic spectra. Besides its remarkable success in explaining the existing experimental 
data in graphene-based systems, the present theoretical approach offers a general scheme for the 
microscopic understanding of Fano-like features in a wide variety of other systems. 



I. INTRODUCTION 

The peculiar properties of single and multilayer 
graphenes make these systems the promising basis for 
the future generation of electronic devices. Within this 
context, the analysis of the spectral properties of the 
phonon anomalies observed by means of different opti- 
cal probes has provided a powerful tool not only for the 
characterization of the samples but also for the investi- 
gation of the underlying scattering mechanisms related 
to the electron-lattice interaction. Large part of the in- 
vestigation along this line has been based on the Ra- 
man spectroscopy* 1 - - — Typical main features under in- 
vestigation within this context were the frequency and 
the linewidth of the phonon anomalies, whose trend as 
a function of doping was found to be in good agree- 
ment with what expected from the theoretical calcula- 
tions of the phonon self-energy, 1 -'! 2 As an alternative 
route, phonon peak anomalies at uj 0.2 eV were de- 
tected also in the mid-infrared optical conductivity of 
bilayer graphenei 13 ' 14 Quite interestingly, unlike in the 
Raman spectroscopy, in this case a strong dependence 
of the phonon peak intensity as well as of its lineshape 
asymmetry on the gate voltage was reported. Under- 
standing and controlling the mechanisms responsible for 
these features at relatively small doping and perpendic- 
ular electric fields is of fundamental importance. 

From a theoretical point of view, the very evidence of a 
strong infrared (IR) phonon activity, as the one reported 
m Refs. [Hill can be considered puzzling in graphenes 
which, having atoms of only one specie (carbon) , present 
a very small intrinsic dipole. In bilayer graphene, for in- 
stance, there are four carbon atoms in the unit cell, as 



depicted in Fig. [T] where atoms Bi and A2 are connected 
by the vertical hopping 71. There are two in-plane op- 
tical modes: an infrared active anti symmetric (A) E u 
mode, which corresponds to (out-of-phase) lattice dis- 
placements in the two layers (Fig. [TJl) ; and a symmetric 
(S) Eg mode, which is associated with in-phase displace- 
ments and is Raman active (Fig. [T|b). If the two layers 
were completely decoupled, all the atoms would be ex- 
actly equivalent to each other. It is thus clear by direct 
inspection that, although allowed by symmetry, the E u 
mode would not induce any dipole and hence it would 
have no IR intensity. Beyond this simple model of decou- 
pled layers, in real systems, the interlayer hopping would 
induce a slight inequivalence between the atoms (Ai, B2) 
and atoms (A2, Bi) and hence a finite electrical dipole 
under the E u lattice displacements. The static dipole as- 
sociated with such physics is however three orders of mag- 
nitude smaller than what is experimentally observed 1 - 3 ", 



E u antisymmetric (A) mode 
A 2 . B 2 



Eg symmetric (S) mode 



A, Bi 



FIG. 1: Atomic structure of bilayer graphene and lattice dis- 
placements for the E u antisymmetric (A) phonon mode and 
for the E g symmetric (S) mode. Labels denote the A and B 
sublattice in each plane. Solid links between atoms represent 
the in-plane 70 hopping, vertical dashed links the interplane 
hopping 71 . 
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FIG. 2: (a) Lowest order contribution to the current-current 
response function, (b) Diagrams involved in the charged- 
phonon effect. Solid and wavy lines are electronic and phonon 
propagators, respectively, white circles are the electron- 
phonon scattering operator while black squares represent here 
the current operator that couples light (dashed lines) to the 
particle-hole excitations. 

so that this effect alone cannot account for the huge in- 
crease of the phonon intensity upon gate-induced doping 
reported in Ref. [TU 

A guideline to understand the origin of this huge 
enhancement comes from the comparison with other 
carbon-based compounds, like fullerenes, where also a 
similar increase of the phonon intensity upon (chemical) 
doping was observed. This effect was explained in those 
materials in terms of a charged-phonon theory 15 i 16 

The basic idea of the charged-phonon model can be 
understood by considering the electronic current (J) re- 
sponse function 

XjM = - J dt(T t j(t)j(0))exp[iut], (1) 

which is related, within the Kubo approach, to the optical 
conductivity as: 



A typical diagram contributing to the electronic back- 
ground at the lowest order is depicted in Fig. [5^, corre- 
sponding to the single-bubble approximation. Effects of 
the electron-phonon interaction on the electronic back- 
ground are commonly taken into account by replacing 
the non-interacting Green's functions in Fig. [5^ with the 
Green's functions evaluated in the presence of electron- 
phonon interaction] 17 ' 18 These processes lead to a smear- 
ing of the optical features and to a redistribution of the 
optical spectral weight. They are however not associated 
with the onset of resonances at the characteristic phonon 
energies in the optical conductivity. 

A most interesting class of diagrams, analyzed by M.J. 
Rice in the context of the charged-phonon effect, is de- 
picted in Fig. [2b- Their contribution to the resulting 
optical conductivity can be described as; 15 ' 16 

&Xjj(u) = IXjphM| 2 DphH, (3) 

where D p h(w) is the phonon propagator of the IR active 
phonon mode considered, and Xj'ph(w) is the "mixed" 
response function between the current and the electron- 
phonon scattering operators (see Sec. [TT] for a more de- 
tailed definition). The resulting optical conductivity of 



this contribution is thus proportional to the phonon prop- 
agator of the optically-coupled lattice vibrations, and it 
presents typical resonances at the corresponding phonon 
frequencies. 

It is worth stressing that the phonon becomes here op- 
tically visible thanks to an intermediate process [xjph(w)] 
where the light couples to particle-hole electronic ex- 
citations. The function Xjph ac ts therefore as a pref- 
actor of the magnitude of the electronically-induced 
phonon resonances. Pristine fullerenes and organic mate- 
rials, for which the charged-phonon theory was originally 
proposed , 15 ' 16 in the absence of doping are semiconduc- 
tors with a band gap significantly larger than the phonon 
energies. In this case the complex function Xjph can be 
reasonably assumed to be real and proportional to the 
charge concentrations n, 

&Xjj(u) oc nD ph (u). (4) 

Eq. Q summarizes in a nut-shell the essence of 
the charged-phonon effect, where the infrared phonon 
activity is triggered-in by the coupling of a lattice 
mode v with the optically-allowed electronic particle-hole 
excitations ! 15 ' 16 As it was shown in Ref. [l9|, the physics 
underlying the charged-phonon effects is intimately re- 
lated to the onset of Fano-like lineshape asymmetries. 20 
Such unified charged-phonon-Fano theory was also em- 
ployed to analyze the spectral properties of infrared op- 
tical phonon in pristine graphite^! and in multilayer 
graphenes with different stacking orders^ 

The purpose of the present paper is to provide a de- 
tailed microscopic derivation of the charged-phonon ef- 
fect in the optical spectroscopy of graphenes. To this 
aim we focus on bilayer graphene as the simplest and 
paradigmatical example. We will show how all the infor- 
mation related to the phonon intensity and Fano asym- 
metry can be evaluated in terms of a unique quantity: the 
current/electron-phonon response function Xjv Within 
this context we evaluate the dependence of the optical 
properties of the phonon resonance on microscopical pa- 
rameters tunable by means of external gating. The cor- 
respondence between the Fano theory and the charged- 
phonon effect is derived microscopically, and a generaliza- 
tion of the charged-phonon effect to the Raman response 
is also provided. 

The structure of the paper is the following. In Sec- 
tion [il] we summarize the main concepts of the charged- 
phonon theory and introduce the mathematical tools em- 
ployed. In Section IIIII we present an analytical discus- 
sion about the correspondence between the Fano and the 
charged-phonon theories, and a suitable quantification of 
the optical properties of the phonon resonances is intro- 
duced. The detailed evaluation of the infrared properties 
in the specific case of ungapped bilayer graphene is pre- 
sented in Section HVl which is generalized in Section [V] 
to the case of gapped bilayer graphene in presence of 
external gate voltage. The role of the electronic struc- 
ture and of the breaking of particle-hole symmetry in 
the band structure is discussed in Section I VII in bilayer 
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graphene and compared with bulk graphite. A gener- 
alization of the charged-phonon theory and of the Fano 
interference analysis for the Raman spectroscopy, within 
an effective mass approximation, is finally presented in 
Section I VIII A summary of the present work and con- 
clusion can be found in Section IVIII1 Appendices [A] and 
iBl provide all the details of the analytical evaluation of 
the charge-phonon theory for IR spectroscopy of the un- 
gapped bilayer graphene in the clean limit, and details 
about a suitable generalization for the Raman response. 



II. CHARGED-PHONON THEORY IN BILAYER 
GRAPHENE 

In order to apply at a quantitative level the concepts 
of the charged-phonon theory in graphenes, in this Sec- 
tion we introduce the electronic band structure and the 
electron-phonon Hamiltonians as well as the relevant re- 
sponse functions which will provide the analytical tools to 
investigate the properties of the phonon peaks in the op- 
tical conductivity and in the Raman response. We focus 
here on the bilayer system as the most simple and repre- 
sentative since the single-layer graphene does not present 
any IR phonon mode. We work in the 4x4 basis of the 
atomic orbitals, as depicted in Fig. [Tk. We introduce the 



four-vector defined as "J 



' I k,<T> ^lk,<r' a 2k<r' ^2k,<r)' 



,t 



where at, and bt operators create an electron with 
spin a in the layer i and on the sublattice A or B, respec- 
tively. 

Considering for simplicity a simple 70-71 model, and 
including a possible asymmetry between the upper and 
lower layer induced by a gate voltage, we can write the 
non-interacting electronic Hamiltonian as: 



Ho 



P,<jHp^p,lT, 



(5) 
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(6) 



70, 71 are the nearest neighbor in-plane and interplane 
tight-binding hopping parameters, respectively, and f p = 
e -i P *a/V3 + 2e lp ^ 2V ^cos(p y a/2). 

Close to the K=(47r/3a, 0) point, writing p = K + k, 
and linearizing around the K point, we can also write: 
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(7) 



Here v = 10 6 m/s is the Fermi velocity for single-layer 
graphene, 71 = 0.39 eV is the interlayer hopping, A is 
the electrostatic energy difference between the layers and 
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FIG. 3: Scheme of the band structure close to the K point. 
The interlayer hopping parameter 71 determines the splitting 
at high energy of the 3 and 4 bands. The interlayer different 
potential A determines the size of the energy gap at the K 
(k = 0) point. 



ir± = h(k x ± iky). The electronic Green's function in 
the imaginary space is thus expressed as a 4 x 4 matrix, 
G(k, iu) n ) = l/[(ihuj n +fi)I — i?k], where /i is the chemical 
potential and hw n — irT(2n+l) are fermionic Matsubara 
frequencies. The electronic bands E^ n are obtained from 
the diagonalization of Eq. (0 , 
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A 2 



J- + — + (hvk) 2 ±\^- + {hvkf W + A^), ( 8 ) 



and are labeled according to Fig. [3] 

We can also define the current operator j a = 

(1/N) J2k a *k, CT ik,a*k,^' where a = x, y, N is the total 
number 01 unit cells, and where 



Jk,, 



e d 
h dk n 



(9) 



In particular, focusing on the current operator along the 
y-axis, from Eqs. ©-([7]) we get 



-evl(a y ), 



(10) 



where A{B) = A®B. 

The electron-phonon Hamiltonian describing electrons 
interacting with the optical modes has also been dis- 
cussed by several authorsi 12 i 23 ' 24 As mentioned above, in 
bilayer graphene there are two optical in-plane phonons 
at q = (see Fig. [T]): the E g symmetric (S) mode and 
the antisymmetric (A) E u one. In the absence of a po- 
tential difference between the two layers, the first one 
is Raman active while the second one is infrared active. 
The relative displacement of the two sublattice atoms A 
and B in the first layer is given for example by: 




r, 2 



■(c q 



iM )e M ( q )e^ r . (11) 



-q,A< 



M is the mass of a carbon atom, uj is the phonon fre- 
quency at the r point, fi — t,l denotes the polarization 
(transverse or longitudinal) and c q , M and ci are the 
phonon creation and annihilation operators respectively. 
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Using q x — (?cos0(q) and q y — gsin^q) the polarization 
vectors e M (q) are given by e;(q) = j(cos </>(q), sin 01 
and e t (q) = i(— sin(/>(q), cos^(q)). Following Ref. 
the interaction between the optical phonons and the elec- 
trons at the K point can be written as 

H int = -V2^ (±1 x u„(r), (12) 

where o- (+) = 1(a), = & x (a), b = 1.42 A is the 

in-plane carbon-carbon nearest neighbor distance, and j3 
is a dimensionless parameter related to the deformation 
potential, whose typical value is (3 — 2.88 i 25 ' 26 In the 
following we shall consider the case of an electric field 
along the y-axis, so that only the lattice vibrations along 
the y direction will couple to the light. Since for q — > 
the result is independent on <fi(q), we take 0(q) = 0, 
so that only the e t polarization vector has a component 
along y. As a consequence we can write, close to the K 
point, the electron-phonon interaction for the v = A,S 
mode as: 



H„ 



where 



E 

k.cr 



* k+q ,^(q)*k^ q - r ( Cq , M + d q)M ), (is) 



V A (q -> 0) = iga z (a x ), 
Vs(q^O) = igl(a x ), 



(14) 
(15) 



and where g = [phv/b 2 ) y/H 2 /2M huj a = 0.27 eV. 
Note that, since e t (q) = — e t (— q), one has that 

liniq^o K(-q) = - lim q ^ K(q). 

We have now all the tools to investigate the full optical 
properties. In this regard it is convenient to define the 
generic propagator for bosonic operators A and B as: 



XAB(iu m ) 



1/T 



cIt(T t A(t)B(0)} exp[iw m T](16) 



where r is the imaginary time, T T the time-ordering oper- 
ator and Hu) m = 2-kttiT the bosonic Matsubara frequency. 

The complex electronic optical conductivity per layer 
a(uS) is obtained from the analytical continuation of 
Xjj(ioJ m ) to the real axis (iui m — >• u> + i0 + ): 



a{u) = - 



y3D 



ifkd 



(17) 



where V 3B = 2dS 2B n with d = 3.35 A being the interlayer 
distance and S 2 ^ = ^/3a 2 /2 the two-dimensional area of 
the graphene unit cell, with a — 2.46 A. 

At the lowest non-interacting order, the current- 
current response function Xjji^m) reduces to the single- 
bubble approximation depicted in Fig. [5Jt: 



Xjj(iu m ) = N s N v —^2Tr\j y G (k 

k,n 



lUJ r 



- iu>m) 



x iyG'o(k,itj„) 



(18) 



where N s — N v — 2 are the spin and valley degeneracies, 
respectively, and Go(k, icj n ) is the 4x4 non-interacting 
electron Green's function which, in the orbital basis, is 
not diagonal. As a consequence, the retarded response 
function Xjji 00 ) wm have contributions coming from both 
the interband and intraband electronic excitations. In 
general it can be decomposed as 



a, (9=1 



(19) 



where a, (3 are the band indexes and describes 
the intraband (a = /?) or interband (a =/= (3) particle-hole 
excitations contributing to the total Xjj ■ The properties 
of the optical conductivity in gated bilayer graphene has 
been theoretically investigated in details in Ref. H?] and 
experimentally confirmed in Refs. [28M3TI One can see 
that, at the phonon energies u ~ 0.2 eV, only the 2-3 
interband transitions contribute to x'jj ( w ) (and therefore 
to the optical electronic background).— Such background 
could be modulated by the charge doping (and therefore 
by the applied gate voltage) so that it drastically vanishes 
for 2\fi\ > 0.2 eV. Within the commonly widespread idea, 
one would expect that this is the electronic background 
that controls the q Fano parameter and then the asymme- 
try of the phonon peak. However, as we shall see below, 
a correct application of Fano theory to graphene leads 
to a different characteristic electronic response function 
related to the Fano effect seen in several experiments. 

Optical properties of the phonon resonances can be in- 
vestigated with the charged-phonon theory by analyzing 
the diagrams in Fig. ^jp. In ungated samples, only one 
phonon, the E u antisymmetric mode, is expected to be 
IR active. We can write thus 



&Xjj{u) = XjA(u)D AA (u)xAu(u), 



(20) 



where XjA is the mixed response function between the 
current operator j y and the electron-phonon scattering 
operator Va- 

However, in the most general case of gated samples 
with A ^ 0, also the symmetric E g mode acquires a 
finite IR activit y 13 ' 14 We can write then 

A X«M = XjA(w)£>AA(w)XAtjM 

+Xjs(w)£>ss(w)x s tj(w) 

+ [x iA (w) J D A sMxstj(w)+h.c.] , (21) 

where v =A, S. 

Eqs. (f2"Tjl can be completed with the Dyson's equation 
for the phonon Green's functions J w 'i 23 ' 24 

[D-\uj)\ vv , = 8 v y [D 1 (uj)] - Xrtv>(u), (22) 

where Dq(lj) = —2frbJo/[h 2 uj 2 — (huj+iO + ) 2 ] is the phonon 
propagator in the absence of electron-phonon interac- 
tion and Xviv'i^) provide the matrix components of the 
phonon self-energy. The quantity luq represents the bare 
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phonon frequencies for a generic mode v in the absence 
of electron-phonon interaction, and it is assumed here to 
be degenerate woa = was = ojq. As it has been discussed 
in detail in Refs. I23ll24]|32l . the mixing between the A and 
S modes (mediated by the self energy xas) is only active 
when A ^ 0. In this case the phonon eigenmodes do not 
correspond any more to symmetric/antisymmetric vibra- 
tions of the atoms in neighboring layers, so that each A 
and S propagator has a double-pole structure, centered 
at the values oj± of the phonon eigenfrequencies. 

Eqs. (j2"T1) - (j2"2")) provide the theoretical tools needed 
to evaluate microscopically the onset and the proper- 
ties of phonon peaks in the optical conductivity. In gen- 
eral, all the information about frequencies and lifetimes 
of phonons is encoded in the phonon self-energy Xv^v' 
whereas the mixed (current / electron-phonon interaction) 
response functions Xjv are related to the intensity and to 
the possible Fano asymmetry of the phonon peaks, as we 
shall discuss in the next Section. 



III. CORRESPONDENCE WITH THE 
STANDARD FANO THEORY 

In order to better clarify the connection between Eq. 
(|2"TT) and the standard Fano theory let us consider first 
the case A = 0, where only the A mode is optically ac- 
tive and the expression (|2Tj) reduces to Eq. (|20[) . In 
addition, for A = the mixed phonon self-energy XAts 
vanishes, which means that the A and S modes coincide 
with the eigenmodes for the lattice vibrations. In partic- 
ular, this implies that the phonon propagator Daa has 
a single resonance, and the main effect of the phonon 
self-energy is to induce a shift of the phonon frequency, 
hu>A — fiujQ + RexAt a( w a) and a finite line broadening 
Ta = — Inr^At a ( w a) ■ Therefore, for oj w oja we can 
approximate the phonon propagator as 



■DaaM 



1 



h(ui - oja) + i^A 



(23) 



Using the relation XAui^) — XjA^), the optical con- 
ductivity can be expressed in terms of the real and imag- 
inary part of the mixed response function XjA and of the 
propagator (f23|) . In particular, introducing the variable 



z = h[u)~LOA}/^ A we have Re£>AA 
and ImD AA = C A A = -1^(1 + 



= D' AA = */r A (i- 

z 2 ), so that 



ImAx 



W 3 a?-W 3 a?~^x' ]A xU 
r A (i + z 2 ) 



(24) 



As a consequence, the real part of the optical conduc- 
tivity (fT7| close to the resonance frequency oja can be 
written asJ^ 



ReAa(w) 



Ia 



q% - 1 + 2q A z 



9i(i + ^ 2 ) 



(25) 



where we defined the Fano parameter q as 



?A 



'x'LM' 



while the prefactor is given by 



[x'ja^a)]- 
oj a T a V 



(26) 



(27) 



As one can see, Eq. (|25j) reproduces the Fano formula, 
where the q parameter controls the asymmetry of the 
peak with respect to a standard Lorentzian profile, that 
is recovered in the limit of q — > oo. As observed al- 
ready in Ref. QjJ the derivation of Eqs. (|2"5|) - ([2l)|) shows 
that the Fano effect stems from a correct implementa- 
tion of the charged-phonon theory. Therefore, for the 
sake of simplicity, in the following we shall regards the 
phonon properties (intensity and lineshape asymmetry) 
arising from this common nature as the "Fano-Rice" ef- 
fect. Moreover, the above set of equations provides a 
general scheme to calculate microscopically the relevant 
parameters that control the shape and the intensity of 
the phonon peak, in particular qA and /a, that are fully 
determined once that the mixed response function XjA is 
computed. 

Before showing explicitly the calculation of XjA we 
would like to make a more direct comparison with the 
standard Fano formalism^ that is often quoted in the 
literature. Following the original work by Fanoj2£ the 
asymmetry parameter q that measures the interference 
effect between a discrete phonon state \<j>) of energy ojo 
and a continuum of electronic states \^)J) can be written 
as 



\T\i) + Pfdoj 
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OJq — OJ 



■u, 0] T\i)nV U0 



(28) 



where P denotes the principal part of the integral, and 
V u measures the hybridization between the phonon and 
the electronic states at the energy w, V u — (ij} u \H\<j>) . 
Here (f\T\i) denotes in general the transition amplitude 
from an initial state |z) and a final state |/). The first 
term in the numerator of Eq. (|28l) represents the response 
of the bare phonon state, i.e. the bare dipole of the sys- 
tem under the lattice distortion, while the second one 
gives the contribution coming from the electronic excita- 
tions. Note that the relevant electronic excitations are 
not restricted to the vicinity of the phonon frequency ujq 
but they involve higher energy states oj' as well. In con- 
trast, the denominator depends solely on the processes 
at oj = ojq and it vanishes if there is no electronic con- 
tinuum at the phonon energy, so that in this case q = oo 
and no Fano asymmetry is expected. In ordinary sys- 
tems the bare phonon intensity {<p\T\i) is large, so that q 
becomes appreciably small only in the presence of a con- 
siderable electron-phonon coupling V Uo , and the observa- 
tion of a pronounced Fano asymmetry is considered as 



(i 



a signature of large electron-phonon interactions £ 3 How- 
ever, the case of graphene is radically different: since here 
the bare phonon activity is negligible, the main phonon 
intensity comes from the particle-hole excitations (Rice 
effect), and the Fano asymmetry can be pronounced even 
in the presence of a relatively small electron-phonon cou- 
pling. In order to show this more explicitly, we introduce 
the complex function: 



where we defined 



did 



, cy) 

u — uj' + irj' 



For (6\T\i) 



0, Eq 



--V u {<ip u \T\i}. 
can be written as: 
X'(^o) 

x"(wo)' 



(29) 



(30) 



(31) 



which reduces to Eq. (|2"oT) derived above from the 
charged-phonon theory. Note that both x' and x" are 
proportional to the electron-phonon matrix element V u 
(i.e. g in our previous notation), so that the strength of 
the electron-phonon interaction cancels out in Eq. (151]) 
as well as in (l26l) . As we shall see below, the response 
function XjA has exactly the form of Eq. ([29|) . and we 
will be able to compute explicitly the function C(tu). 

Few more final observations are in order concerning 
Eqs. ([25]) - ([33]) . First, in contrast to the ordinary case 
where also the bare dipole charge of the phonon must be 
considered, in the case of graphene the two quantities qx, 
I A are not independent. This is again a consequence of 
a correct implementation of the charged-phonon effect, 
i.e. of the fact that in graphene the optical activity of 
the phonon is fully borrowed from the electronic excita- 
tions. This also permits to quantify the "strength" of a 
phonon resonance in the optical conductivity on a more 
rigorous ground. To this aim, for a given phonon mode 
v, two typical quantities are considered in literature: the 
integrated spectral area 



W' v = J dujReAa(uj), 



(32) 



and the phonon intensity W v , as obtained from the Fano- 
like fit in Eq. (gSJ) as 



W v = 7rr„J„. 



(33) 



The quantity W v is considered to be the "bare" spectral 
intensity that the mode would have in the absence of 
the Fano interference. It is straightforward to show that 
these two quantities are related by the formula 



wL = i 



i 



W v , 



(34) 



so that W — > W v when \q v \ — > oo for a symmetric 
Lorentzian profile (Fig. [JJi)- 
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FIG. 4: Sketch of the optical properties of a phonon reso- 
nance in different characteristic regimes: (a) positive sym- 
metric Lorentzian peak (\q\ — oo; |x"ji/(Wi>)| = 0), where 
W = W'; (b) weak asymmetric Fano profile (ex.: \q\ — 5; 
\x'jv(wv)\ = 5|x"jv(w„)|); (c) highly asymmetric case (\q\ = 1; 
1x^(^)1 — Ix"3v( w ^)I)j where the integrated area W' = 
because of the cancellation of positive and negative spectral 
regions; (d) negative phonon peak (\q\ — 0; x ~ 0), where 
the "bare" intensity vanishes, W = 0, although a phonon 
anomaly is visible in the imaginary part of x- The total 
phonon strength P is the same in all the cases. 



Both W v and W' v can be expressed in the common 
Fano-charged-phonon framework via the mixed response 
function Xjv- 



Wv = 



0J V 



and 



wi 



[x?>„)-x£M] 



(35) 



(36) 



It should be however stressed that none of these two 
quantities, W , W, can provide a satisfactory quantifi- 
cation of the magnitude of the phonon optical anomaly. 
For instance, in the case \q u \ ~ 1, where the phonon 
peak asymmetry is strongest and which corresponds to 
|X^v( w f)l = Ix'jviUv)], we get W' v = due to the cancel- 
lation of positive and negative spectral regions [Fig. 
On the other hand, when Xjvi^v) ~ 0, the spectral prop- 
erties are characterized by a sizable negative Lorentzian 
peak whose intensity is driven by x"„(w„) whereas the 
estimated "bare" intensity is vanishing small, W u ~ 
[Fig. EH]. 

Both definitions W and W thus fail to describe the 
actual magnitude of the optical phonon resonance inde- 
pendently of its Fano-like properties. This problem can 
be solved however thanks to the microscopical identifi- 
cation of the optical properties in terms of the mixed 
response function Xjv The simple identification of W 



7 



and W as W v oc x%, W{, oc x% - x"l suggests us to 
introduce a strictly positively defined quantity, 



and where 



(37) 



which we refer to as phonon strength. Note that P — > 
only when both x'jviX'jv ~~ * 0, i.e. when the phonon 
features are indeed vanishingly small. In the following, 
when discussing the magnitude of a optical phonon reso- 
nance, we shall therefore refer to this quantity P v , which 
permits to describe the actual visibility of a phonon struc- 
ture independently of its Fano shape. The robustness of 
the parameter P v to characterize the magnitude of the 
phonon resonance independently of its lineshape and its 
Fano properties is demnostrated in Fig. 0] where all the 
panel where evaluated for fixed phonon strength. 

As a last point we would like to stress that it is the 
mixed response function XjA, and not the electronic back- 
ground Xjji as ^ was considered in Refs. Il5lll6l that de- 
termines the phonon strength and the asymmetry. As we 
shall see this makes a crucial difference. Indeed, while all 
the possible electronic excitations contribute to Xjj 1 on ly 
a subset of them enters into the mixed response function 
XjA (and XjS m the case A / 0), determining in this 
way the exact selection rules for the phonon activity in 
bilayer graphene. 



IV. FANO-RICE PROPERTIES IN UNGAPPED 
BILAYER GRAPHENE (A = 0) 

In the case of unbiased graphene the mixed response 
function XjA can be computed analytically in the bare- 
bubble approximation, namely: 



xV v G(k,iuj n ) , 



iui n + iui m ) 



(38) 



where we use in Eq. (|38l) the bare electronic Green's func- 
tions. Some additional details of the calculation are given 
in Appendix |XJ Due to the multiband structure of the 
system, XjA has the typical structure of a particle-hole 
Lindhard response function, with proper coherence fac- 
tors C™^™ weighting the contributions of the various exci- 
tations between the n and m bands. In particular, using 
the explicit matrix expressions of the j and Va operators, 
one gets 

XjaM = x]aM + x]aH - XmM - X, 3 lH-(39) 

Here n, m are the band indexes and Xj'Ti 10 ) = 7r j l ™( w ) — 
71-^4™ (w), where 

nm/,,1 1 f( E W,n - Mj - f (E^rn - Mj , n s 



c 



JAM 



gevN s N v 



7i 



(41) 



for (n, m) = (1,2), (1,3), (2,4), (3.4) and zero otherwise. 
In Eq. (|4T)j) fi represents the chemical potential, and 77 
is a phenomenological parameter that accounts for the 
damping in the electronic states, so that the clean limit 
corresponds to ij — > OM- Note that, once we identify 



m,n k 

x[f(E Kn - v)- f(E Km ~ fi)}, (42) 



the function XjA(w) can be written exactly as in Eq. 
(|29[) above, with the particle-hole energy difference u/ = 
Ek, n — -Ek.m playing the role of the electronic continuum 
in the Fano theory. Thus, our approach allows one to 
identify the optical intensity of the phonon peaks and 
their Fano parameters in terms of the real and imaginary 
parts of a specific response function, that can be cal- 
culated using the standard diagrammatic theory. This 
means in particular that we can determine: (i) which 
phonon is coupled to the current, and (ii) which elec- 
tronic excitations couple to each phonon mode. 

All these elements can be quantified in an analyti- 
cal way for the A phonon in the case of no gap. First 
of all, note that for fi — the term x]a cancels pre- 
cisely with XjA' whereas xfk^ X^a are both vanishing 
because particle-hole excitations between completely full 
or empty bands are not allowed. As a consequence 
XjAiy) = 0, implying that the intensity of the A mode, 
although it has the correct symmetry, vanishes in un- 
doped bilayer graphene.— As we shall see, this property 
holds true even in the presence of an interlayer difference 
potential A. 

The total mixed response function XjA(^>) as well as 
its single contributions x"^ n (u;) can be easily evaluated in 
the clean limit, r\ = 0, and in the linear approximation for 
the band dispersion, corresponding to the Hamiltonian 
given by Eq. ([7]). Details of the computation are also 
reported in Appendix El It is convenient to introduce 
the dimensionless quantity £(w): 

where A = ge^S™jA-Kh 2 v. At T — we obtain the 
analytical expressions: 



in 



(7i + fiw)(7i-^ + 2|/-i|) 



(7i 

0(|M|-7l)ln 
4MlM|0(7i " 



fiw)( 7 i + hw + 2\n\)_ 
(71 -I- fa;) (71 - 



hio - 2\fx\) 



(71 — fat;) (71 — hjj - 

ImI) + 0(Im|-7i)] 



2| M | 



(43) 



( 7l 2 - ft 2 ^ 2 ) 

7r{0(fiM- 7l )0(/iM-2H+7i) 
-0(h\u\ - 2| M | - 7l ) 

+2[H0( 7 i - ImI) + 9Qn\ - 7i)W%l - nm. 
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FIG. 6: Dimensionless quantities £'(oja), £,"{^a), as well as 
the spectral weights Wa, W' a , Pa and the Fano factor qa as 
functions of the charge concentration n for the clean limit 
(rj — 0) and for finite damping 77 = 5, 10, 15, 20 meV. Curves 
for £'(wa), Wa, W' a , Pa at zero and finite 77 are barely dis- 
tinguishable. 



The real and imaginary part of Xj'a(w) for a representa- 
tive case n = 5x 10 12 cm" 2 (/i = 0.13 eV) are shown in 
Fig. [5k., whereas in panels (b)-(d) the contribution of the 
single interband transitions is also shown. It is interest- 
ing to remark that for |w| > 71 + 2|/j| the term Imx] A (w) 
cancels exactly with lmx 2 A {u), so that only a limited 
low- frequency energy- window 71 < |wj < 71 + 2|/x| con- 
tributes to the imaginary part (and hence, via Kramers- 
Kronig relations, to the real part) of the mixed response 
function XjA ■ This fact guarantees that the result is valid 
even if the calculation has been done using the linearized 
graphene bands instead of the full periodic band struc- 
ture. 

According to Eq. (|2~5"1) . the computation of the opti- 
cal phonon spectra requires the evaluation of the quan- 
tities P3")) - (fP| at the phonon frequency wa- In order to 
determine wa, one should in principle solve the phonon 
Dyson's equation (f2"2"j). However, the shift of the renor- 
malized phonon frequencies u) v is only a few meVs, which 
is much smaller than wo 200 meV, so that we can re- 
place wa with wo- In any case, as one can see in Eq. (|39[) . 
since the low-energy 2-3 interband transitions are not al- 
lowed in XjA, the lowest threshold for the particle-hole 
excitations is determined in the clean limit 77 = by the 
1-3 interband transitions with uj > 71 = 0.39 eV, so that 
ImXjA^A) = 0. According to Eq. (|2"o]l this implies that 
the A phonon peak in the clean limit has no Fano asym- 
metry (I^aI = 00), and the bare phonon peak intensity 
Wa (|33| coincides with the integrated area (Wa = W' A ) 



[see Eq.flSU)]. From Eq. Q we get : 

= ^5d£'Va) = Acto^'Va), (45) 

where cto = e 2 /4hd « 1816 fi" 1 cm" 1 and A = 
(^/3/ir)g 2 /(Hv/a) 2 is the dimensionless phonon coupling 
which, using the value g = 0.27 eV estimated in Section 
HI1 results in A = 6 x lO" 3 ^ in agreement with the exper- 
imental estimates given in Refs.l5Hl3l. Finally, we obtain 
that nA 2 /oj A V 3D = 4.35 x 10 3 ft^cm" 2 . 

To elucidate the doping dependence of the optical 
properties of the A mode, we show in Fig. [6] the quan- 
tities £'(wa), £"(wa) as functions of the charge doping 
n, along with the parameters Wa, W' a , Pa and qa- As 
mentioned above, due to the absence of 2-3 interband 
transitions, the imaginary part £(w) is zero at w = co A 
in the clean limit 77 = and relatively small for finite 77 
(Fig. [6p), so that the spectral weight is dominated by 
the real part £(uj a ). In this context the small contribu- 
tion at u — lja for finite 77 coming from the gapped 2-4 
and 3-4 interband excitations gives rise to a finite Fano 
factor |qa| < 00 (Fig. |6f) but it does not affect sensibly 
the spectral weights Wa, W' a and the spectral strength 
Pa- In particular, the doping dependence of these latter 
quantities, in the range n here considered, is dominated 
by the 3-4 interband transitions [the third term in Eq. 
(|35])], where £'(w A ) oc so that, for £"(w A ) < £'(w A ), 
W A = W' A = P A oc K>a)] 2 <x | m | 2 « n. 

Note also that, although the imaginary part £"(wa) 
does not contribute to the spectral weights Wa, W' a and 
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Pa, it determines the magnitude of the Fano asymmetry 
factor, as shown in Fig. [6f. As mentioned above, this 
is triggered by the finite spectral weight in £"(cja) due 
to the broadening r\ of the higher-energy transitions. In- 
creasing the charge concentration n, and hence the chem- 
ical potential leads to an overall increase of £"(u>a) 
and, since £'(wa) is related by the Kramers-Kronig rela- 
tions to the low energy part of to a similar increase 
of £'(wa)- The charge-doping dependence of the magni- 
tude of £'(cl>a) an d f (wa) is thus similar, making the 
Fano factor (/a almost independent of n. As we are go- 
ing to see in the next Section, the situation is different 
for the symmetric S mode, when its infrared activity is 
triggered by a finite energy difference A^O. 



V. FANO-RICE PROPERTIES IN GAPPED 
(A / 0) BILAYER GRAPHENE 

In the above Section we have addressed the opti- 
cal properties of the antisymmetric E u mode in bilayer 
graphene in the absence of any electrostatic potential 
gradient A between the two layers. However, in most 
cases, the gating of the samples as well as the influence 
of the substrate give rise to a finite potential difference 
between the two layers. In this case, the antisymmet- 
ric A and symmetric S modes depicted in Fig. [1] are no 
more eigenvectors of the lattice dynamics although they 
still represent a suitable basis to investigate the optical 
properties of bilayer graphene. According to Eq. ([21]) we 
can distinguish three different channels responsible for 
the onset of phonon peaks in the infrared conductivity: 
the direct coupling of the electron current with the A 
mode, which is already present for A = [first term of 
Eq. (j2"Tj) ]; the direct coupling with the S mode, which is 
induced by the presence of the gap A and vanishes for 
A — » [second term of Eq. ([2T])]: the mixed A-S mode 
optical coupling where the incoming light first excites the 
A lattice vibrations, which develop a S component due 
to the hybridized A-S phonon self-energy, and finally the 
light is re-emitted through the coupling of the A mode 
to the current [third term of Eq. (|2~Tj)]. 

A. Undoped case n = 0, A / 

An interesting case to elucidate the role of the gap in 
triggering on the S-channel is the undoped case (fi = 
0) where the antisymmetric mode E u is not involved 
[x"a( w ) = 0]r^ and the only coupled lattice mode is 
the symmetric E g one induced by the finite A. Such 
condition has been experimentally realized in Ref. Il4l . 
In Fig. [7^, we show, for different values of A, the 
real and imaginary part of the dimensionless quantity 
£s( w ) = Xis(w)/A, which rules the optical properties of 
the phonon resonance. An important difference here is 
that the low-energy interband transitions 2-3, that were 
forbidden in XjA, are allowed in the response function 
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FIG. 7: (a) Real part (upper panel) and imaginary part (lower 
panel) of the dimensionless response functions £(oj) for fi — 
and for different values of the gap: A = 0.05, 0.1, 0.2 eV. Here 
we set T = 10 K and r\ — 10 meV. (b) Gapped dispersion (left 
panel) and corresponding density of states N(e) (right panel) 
for A = 0.2 eV. The vertical arrows mark the gap A at k = 
(K point) and the optical gap A at the bottom of the mexican 
hat. In this case the DOS presents a divergent behavior. 

XjS. The imaginary part of £s(w) is thus finite at low en- 
ergies (limited only by the opening of the optical gap A) 
providing a finite Fano interference as long as ws < A. 
The sharp features in £g (w) are reflected in peaked struc- 
tures in the real part £g(w). In particular, the structure 
at u) = 0.4 — 0.5 eV can be associated with the transitions 
between bands 2 and 4, with lowest characteristic energy 
at w 7i + A/2, whereas the structure at low energies 
oj = 0.05 — 0.2 eV reflects the opening of the band gap on 
the transitions 2-3. It is worth to note that, due to the 
mexican-hat shape of the electronic dispersion, the lowest 
energy threshold is not determined by the gap A at the K 
point but by the actual optical gap A = A71 / \J~i\ + A 2 
that lies at a finite momentum k at the bottom of the 
mexican hat (Fig. [7}d). The reduced dimensionality of 
the electronic dispersion in these points gives rise to a 
singular behavior in the density of states, which is re- 
flected in a corresponding behavior in the particle-hole 
excitations in £g(a;). Such peaked structure in £g(w) has 
important consequences also on the real part of £s( w )i 
resulting in a strong peak at uj w A. When A w 0.2 
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eV (A « 0.233 eV) both ^(w) and ^'(w) are peaked 
at u) rj cug and the phonon strength is expected to be 
strongly enhanced. 

This trend, as well as the dependence on A of all the 
optical properties of the S phonon peak, are shown in 



Figs. E For A < 



LOS 



0.2 eV the phonon strength 



Ps is increasing with A signalizing the switch-on of the 
symmetric E g mode. Such increase of the strength is 
mainly driven by £"(ws) while £'(ws) ~ 0. This results in 
a negative peak where Ws ~ and a negative integrate 
area W$ < 0. The vanishing of £' (cos) results also in 
a Fano factor qs ~ 0, which indeed corresponds to a 
negative symmetric shape. 



B. 



Doped case n ^ 0, A^O 



Unless using a double-gate device^ where the top and 
bottom gates are tuned to set the doping or the gap to 
zero, in single gated devices the gating induces at the 
same time an inversion symmetry breaking and a finite 
doping. In this most common case, both the E u and E g 
modes are simultaneously IR active. In order to inves- 
tigate theoretically the optical properties of the phonon 
resonances, one has to employ the full Eq. ([2Tj). where 
both modes are presents. The phonon spectral proper- 
ties are then much more complex than in the A = case. 
Not only we have different phonon channels contributing 
to the total features, but also the phonon propagator of 
each channel [e.g. the Pss(w) propagator], under partic- 
ular conditions, can develop a double-pole structure, as 
discussed in Refs. 5.23.24.32 in the context of the Raman 



> 




AS 



FIG. 9: Phase diagram of the relative relevance of different 
phonon channels _Ras = (Pa — Ps)/(Pa + Ps)- In the green 
region optical features associated with the antisymmetric E u 
mode are dominant, in the red one the dominant mode is the 
symmetric E g one. Also shown here is the location in the /i-A 
space of the experimental available data [Refs. [3 (squares), 
[hi (circles)]. Adapted from Ref. fl9l . 



spectroscopy. 

In general, we can attribute a different role to the sev- 
eral quantities appearing in Eq. (|21[) . Here the mixed 
response functions Xjv are mainly responsible for the 
magnitude and the Fano lineshape of the phonon fea- 
tures, while the pole structure of the phonon propa- 
gators £>„„', is associated with the frequencies of the 
phonon resonances and to their linewidth. Keeping in 
mind this distinction, and since the frequency structure 
of Xjv(w) varies over electronic energy scales while the 
phonon self-energy gives rise to a splitting of the S and 
A mode frequencies of few cm -1 , we can in a very good 
approximation evaluate the functions Xjufa) a ^ the bare 
phonon frequency u) ~ loq. In this context, the rele- 
vance of each phonon mode in the infrared spectroscopy 
is ruled by the characteristic phonon strength of the cor- 
responding channel, i.e. Pa = 7r|XiA(wo)| 2 /wo 1 ^ for the 
antisymmetric E u mode, associated with the first term 
in the r.h.s. of Eq. (|2T|). and P s = n\Xis{uo)\ 2 / uqV 
for the symmetric E g mode, associated with the second 
term in the r.h.s. of Eq. (|2"Tj) . It is also possible to de- 
fine a mixed channel, related with the third term in the 
r.h.s. of Eq. (|2i"j) , characterized by a phonon strength 
Pas = 7r VI^' A ( CJ °)^sti(wo)|/woU. This channel is how- 
ever quite weak in most cases ^ and we neglect it in the 
following discussion. 

The comparison of the phonon strengths associated 
with different A and S phonon modes is shown in Fig. 
[9j where we plot the relative intensity Pas = (Pa — 
Ps)/(Pa + Ps) as a function of the chemical potential fi 
(doping) and the band gap A induced by the vertical elec- 
tric field. As we discussed above, Pa is essentially driven 
by the doping, whereas Ps is induced by the gap A. Ac- 
cording to the relative position in the /i vs. A phase 
diagram, we can thus predict one mode to be dominant 
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with respect to the other one. A "phonon switching", 
namely the switch from the dominance of one phonon 
mode to another mode, is possible as a function of the 
gate voltage. In Fig. [9] we also show the location in 
the fjL-A space of the experimental available data from 
Refs. [Ulli While the data of Ref. OJ, in the double- 
gated device, were mainly collected along the neutral line 
/i = 0, the optical conductivity measurements of Ref. [l3| 
in the single-gate geometry span a much wider region, 
going from regions where the antisymmetric mode E u is 
expected to dominate to regions where the optical fea- 
tures can be attributed mainly to the symmetric mode 
Eg, The actual evidence of such theoretically-predicted 
phonon switching in the experimental data was discussed 
in Ref. [HI, where we refer the reader for more details. 

One should underline that the phonon switching is not 
directly related to the appearing of a double-peak struc- 
ture in the optical conductivity, as observed in Raman 
spectroscop y 5 ' 23 ' 24 ! 32 What we are describing here is the 
dominant intensity of one phonon with respect to the 
other one, in other words, which phonon mode is most 
coupled to the light. This information is encoded in the 
mixed response function Xjv which describes the coupling 
of the light with the mode v. In this respect, as shown 
in Fig. [3J in bilayer graphene the two channels are to 
a good extent mutually exclusive. A different matter is 
the possibility, within a given phonon channel, to develop 
a two-peak phenomenology. This issue is related to the 
presence of a large off-diagonal phonon self-energy xa.s of 
the same order of the diagonal phonon self-energies xa,a, 
Xs,s in the 2x2 space of the phonon modes v =A, S. On 
this respect, since the phonon self-energy is a quantity 
which is shared in both Raman and IR spectroscopy, we 
expect that phonon features in optical conductivity can 
develop double-peak structures in the same regions where 
the Raman spectroscopy sees them. Things are however 
more complicated in the optical conductivity case since, 
as we show in Section IVIIl Raman spectroscopy in bi- 
layer graphene is dominated only by the S channel. On 
the other hand, in the optical conductivity, we expect 
that the condition xa,s 7^ 0, ruling the double-peak fea- 
tures, would arise in the same phase space where the 
mixed channel Pas [second line in Eq. (|2"Tj) ] is of the 
order of Pa, Ps- I n this situation spectral interferences 
between the different channels can occur making the sce- 
nario more complex than in Raman spectroscopy. 



VI. FULL TIGHT-BINDING MODEL 

The simplified model considering only the leading tight 
binding terms 70 and 71 has some significant limitations. 
For example, within this model the intensity of the IR 
active E u drops exactly to zero in the undoped limit 
/i — > 0. The generalization of this model to bulk graphite 
would predict thus no IR activity at all, despite the clear 
evidence of a phonon resonance reported already in the 
70s i 36 ' 37 As a matter of fact, the evidence of such phonon 



intensity in graphite was reconciled with the charged- 
phonon theory in Ref. [2l] where it was shown that the 
inclusion of higher order tight-binding terms, in partic- 
ular of the ones breaking the particle-hole symmetry, is 
responsible for the observed phonon activity. It is thus 
interesting to investigate to which extent these higher- 
order tight-binding terms can affect the results for the 
bilayer graphene. 

We address this issue by including explicitly in the 
Hamiltonian of bilayer graphene the higher-order tight- 
binding hoppings 73 , 74 as well as the crystal field 5 which 
differentiate the atoms Bl, A2 from the atoms Al, B2. 
Close to the K point, we can write: 
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«47T_ 


V 3 TT + 


V 7T + 


<5 + A/2 


71 


V4ir_ 


V4TT+ 


7i 


S- A/2 


VTT- 


V «37T- 


v 4 ir + 


vir + 


-A/2 



(46) 



where = 71/70 ■ We consider here typical values of 
73 = 0.29 eV, 74 = -0.13 eV and S = 0.022 eV^ We 
evaluate consequently also the current operator 



e d 
h dk„, 



= -evl{(jy) - ev 4 a y (I) + ev 3 — yy (4 7) 

and the electron-phonon scattering matrices 

V A = iga z {vx) + i940- x (a z ), (48) 

V s = igl{a x )-ig 3 g ' . (49) 

The terms g 3l g^ represent the electron-phonon coupling 
associated respectively with the hopping terms 73 and 74 
and they can be related to the corresponding deformation 
potentials which have been recently evaluated by DFT 
calculations* 3 ^ We get namely g 3 = 0.033 eV and 34 = 
0.018 eV. 

As two representative limits, we show in Fig. ITUa.b 
the phonon strength Pa and the Fano parameter gA for 
the E u mode in the ungapped case as functions of the 
doping n, and in Fig. UOb .d the phonon strength Pg and 
the Fano parameter gg for the E g mode in the undoped 
case as functions of the electric field parametrized by the 
gap A at the K point. As found for bulk graphite in 
Ref. [U a finite phonon strength P A ~ 0.62 x 10 3 Q.' 1 
cm~ 2 (weakly visible on the scale of Fig. ITUa ) is now 
triggered by the higher-order tight-binding terms at the 
neutral point. However, apart from a slight asymme- 
try for electron and hole doping, the inclusion of such 
tight-binding terms does not change qualitatively the re- 
sults at finite doping, with a roughly linear increasing 
of the phonon strength Pa as function of \n\. In addi- 
tion, such residual phonon strength P\ at the neutral- 
ity point is easily overwhelmed in gated systems by the 
presence of the symmetric E g mode with Pg » Pa (Fig. 
[TUb). The A dependence of the optical properties of such 
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FIG. 10: Phonon strength Pa (panel a) and Fano parameter 
q\ (panel b) for the E u mode as functions of doping n for 
A = including all the relevant tight-binding parameters 
(black solid line) and using the simple 70-71 model of Section 
4 (red dashed line), with T — 10 K and r\ = 20 meV. Panels 
c and d: same quantities Ps, <Zs for the E g mode at n = as 
functions of the gap A. 



mode (phonon strength P§, Fano parameter qg) is also 
barely affected by the presence of the higher-order tight- 
binding terms, so that we can conclude that the scenario 
presented in Sections IIV1 IVl is generally robust in realis- 
tic materials against the inclusion of higher order tight- 
binding terms, except in the case of ungated and undoped 
bilayer graphene, as it could be relevant in suspended 
samples. 



VII. CHARGED PHONON THEORY FOR THE 
RAMAN RESPONSE 

In the previous Sections we have outlined in de- 
tails how a quantitative implementation of the charged- 
phonon theory gives rise in the infrared optical spec- 
troscopy of bilayer graphene to a phonon intensity and 
to a asymmetric Fano-like shape for the A mode as well 
as for the S mode once this latter is turned on by the 
interlayer potential. In particular, such analysis shows 
that for gate potential V g close to the charged neutrality 
point, the negative phonon peak can been attributed to 
the S mode with a phonon activity strongly dependent 
on the gate voltage V g and with a strong Fano charac- 
ter with qs 0, induced by the interference with the 
low-energy 2-3 interband particle-hole excitations. 

One can wonder why the same S mode, under the 
same conditions (namely tuning charge concentration 
and bandgap A) does not present in the Raman response 
any significant Fano asymmetry and strong dependence 
of the phonon intensity upon V g &Z& In this Section we 
show that this different behavior for the infrared and Ra- 
man spectroscopy can be naturally explained within the 
context of the charged-phonon theory. 

In our analysis we focus on the possible changes in 



the Raman features of the phonon resonances at ~ 1590 
cm -1 as functions the charge doping and the bandgap, 
i.e. quantities that typically affect only the low-energy 
excitations close to the Fermi level. We do not address 
the possible dependence of the phonon intensity on the 
external laser energy 40 ! 41 We also assume a relatively 
small gate-induced doping, so that the chemical potential 
is much less than the laser energy, in a region where the 
absolute Raman intensity of the phonon peak at ~ 1590 
cm -1 is constant^ Within this context, we employ the 
effective mass approximation^ when the only relevant 
quantity is the Raman shift u> = uJi n — i^mit, of the order 
of the phonon energy, u> ~ ujq- Refs. I40il4ll showed also 
how the main electronic transitions responsible for the 
phonon intensity were associated with high energy pro- 
cesses, of the order of the laser energy uj- ln ~ 1 eV or of 
the order of the 7r-bandwidth W ~ 6 — 7 eV, in any case 
much larger than the phonon energy scale. As we shall 
see, we will recover this result in our simplified scheme 
and we will show how it can explain the different phe- 
nomenology of Raman spectroscopy with respect to the 
infrared spectroscopy. For simplicity we consider unpo- 
larized isotropic Raman scattering. The Raman intensity 
I(uj) can be related to the imaginary part of the Raman 
response function^ namely 

J(w) = --[I + 6(w/T)]ImxRM(w), (50) 

where b(x) — l/[exp(x) — 1] is the Bose-Einstein function 
and xrm(w) is the analytical continuation to the real fre- 
quency axis of the Raman response function XRM(*w m ), 
whose explicit definition and derivation are given in Ap- 
pendix [B] Just as in the case of the optical conductivity, 
the irreducible part of xrm(^) provides the electronic 
Raman background, which is proportional to the current- 
current response function involved in the optical conduc- 
tivity. In a similar way, within the charged-phonon the- 
ory, the bubble diagrams mediated by a phonon propa- 
gator, as in Fig. O3, are associated with the onset of the 
phonon peaks in the Raman response. In particular, we 
can write, for A^0, 

Axrm(w) = x 7 sM-Dss(w)xs+ 7 H 

+X 7 a H £>a a H x a 1 7 M 

+ [X7sH-DsAHxAt 7 H 



h.c] ,(51) 

where X7i/( w ) are Raman mixed response functions in- 
volving one Raman vertex 7, described by the effective 
mass approximation, and one electron-phonon scattering 
operator. Explicit expressions for the Raman vertices 
for different polarizations are given in Appendix [B] It 
is worthwhile to point out that close to the K point the 
Raman vertex scattering operator for the xx polarization 
involved in Eq. (ISTl) reads 



(52) 



7£ x oc I(a x ) oc iV s . 



This observation permits us to relate the Raman spec- 
troscopy in the effective-mass approximation to the re- 
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sponse functions previously discussed. In particular, us- 
ing the relation (|52|). we can relate the Raman mixed 
response functions to the phonon self-energy. We obtain 
for instance x 7 sM oc xssM and x 7 aM oc xsa(w). 
Just as in the infrared response, also in this case the Ra- 
man spectral strength of each mode will be pR aman oc 
|x 7 i/(w l/ )| 2 and the Fano asymmetry factor gR aman = 

A crucial difference with respect to the infrared spec- 
troscopy is the different dependence of the Raman mixed 
response functions Xiv on the high-energy interband 
transitions when compared to the infrared response func- 
tions Xju ■ One should keep in mind that the linear Dirac- 
like dispersion e = hv\k\ of the single-layer graphene ex- 
tends only up to 3-5 eVs. However, as we saw in Sec- 
tion IIV| in the clean limit and in the absence of the 
gap A the imaginary part of x.ja(w), and hence also 
the real part due to the Kramers-Kronig relations, are 
uniquely determined by low-energy interband transitions 
7 < M < 7 + 2|/i|, so that the high-energy cut-off E c 
does not play any role. One can see that such results hold 
true also for Xjs(u) and for generic finite i] and A, in a 
sense that that the imaginary part of the infrared mixed 
response functions Xjvi*^) is not divergent for lo — > oo 
so that the high-energy cut-off E c can be safely set to 
infinity. 

The situation is different when the Raman response 
function is considered. In this case, using Eq. (|52[) . 
we get x 7 S oc xss and x 7 A oc xsa- The function 
Xss(w), which corresponds to the phonon self-energy for 
the symmetric E g mode, has been widely analyzed in 
the literature . 12 ' 23 ' 24 In particular, its imaginary part in- 
creases linearly at high energy, Xss( w ) ~ u - As a con- 
sequence, from the Kramers-Kronig relations, the mag- 
nitude of the real part of xss(w) is dominated by the 
high-energy processes, i.e. by high-energy cut-off E c , 
Xgg(w) ~ E c . On a physical ground, this cut-off E c can 
be identified with the 7r-bandwidth W or with the highest 
energy uj- m above which the effective mass approximation 
breaks down. 

From a careful inspection, one can see that both the 
real and imaginary parts of almost all the mixed re- 
sponse functions x-y,v are only weakly dependent on the 
high energy cut-off E c . The only exception is the real 
part of the Raman coupling with the S mode, Rex 7 ,s, 
which scales with E c . Since this is the highest energy 
scale involved in the system, we have x 7 s ^ X 7 S' X 7 A' 
X 7 a- A s a consequence we see that the phonon Raman 
spectroscopy is dominated by direct coupling of the Ra- 
man scattering operator with the E g symmetric mode 
[first line of Eq. (|5"T|) ]. whereas other channels involv- 
ing the E u are marginal, independently of the charge 
doping level or of the interlayer potential difference A. 
Other interesting consequences: i) since x 7 s ^ X 7 s> tne 
Fano factor for the dominant S channel results |gR aman | — 
|x 7 s( w s)/x 7 s( a; s)| ^ 1 and the Raman phonon peaks 
are expected to be positive and symmetric, in agreement 
with the experiments; ii) the magnitude of the Raman 



phonon strength pR aman & [x' jS ] 2 is mainly dominated 
by the energy cut-off E c , so that it does not depend sig- 
nificantly on the low-energy features related to the charge 
doping or to the opening of the gap A, also in agreement 
with the experimental observations. 

Note that the fact that the phonon Raman response 
is dominated by the direct coupling of polarizability to 
the S mode does not exclude that a double-peak structure 
could be observed. Indeed, as we saw in Section|Vl under 
suitable conditions the phonon propagator Z?ss (w) itself, 
as well as £>aa(w), can develop a double-peak structure 
as a consequence of the phonon-modes hybridization trig- 
gered by the mixed phonon self-energy XAS^ ' 2ai24 ' 32 



VIII. CONCLUSIONS 

In this paper we have provided a comprehensive deriva- 
tion of the charged-phonon theory applied to investigate 
the optical properties of the phonon peaks in the opti- 
cal conductivity of bilayer graphene. The origin of the 
phonon activity and its relation with the occurrence of 
a Fano effect have been elucidated. The dependence of 
these properties on the tunable microscopical parame- 
ters, i.e. the doping and the bandgap induced by an 
external gate voltage, has been discussed in detail. We 
have also compared the charged-phonon theory in the in- 
frared and Raman spectroscopy, accounting for the dif- 
ferent phenomenology of the phonon peaks in these two 
different optical probes. 

The theory presented here provides a suitable tool to 
characterize quantitatively bilayer graphene in terms of 
the intensity and Fano asymmetry of the infrared phonon 
peaks. Further future developments of the present analy- 
sis could investigate the dependence of the phonon opti- 
cal properties on the different symmetry breakings* 4 ^ - — 
in order to provide a fingerprint for the possible un- 
derlying instabilities. The microscopical characteriza- 
tion of the phonon optical properties sheds also light 
on the underlying physics of the electron-phonon inter- 
action. While the present work was focused on bilayer 
graphene, as the simplest graphitic system with infrared 
active modes, the theory presented here can be general- 
ized in a straightforward way to multilayer graphene, 22 
and to bulk graphite^ 

The general scheme discussed in this work to investi- 
gate within a microscopic many-body approach the Fano 
effect can be applied, with the due modifications, also 
to different classes of materials. One remarkable exam- 
ple is provided by layered systems with different atomic 
species in the units cell, for instance M0S2, Bi 2 Se3, where 
phonon anomalies have been detected in infrared and 
Raman spectroscopy 5 1 ' 52 A second interesting case is 
the As phonon mode in the pnictide Ba(Fei_ a; Co^ AS2, 
which displays a concomitant presence of intensity vari- 
ation and asymmetry as a function of doping and/or 
temperature, in particular across the magnetostructural 
transition. 53 In all these cases the understanding of the 
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phonon anomalies can shed new light on the underly- 
ing bulk electronic structure. In this sense phonon spec- 
troscopy can represent a powerful and alternative tool to 
investigate electronic excitations also in correlated mate- 
rials, provided that a correct implementation of the Fano- 
Rice theory is used to relate the phonon and electronic 
response, taking into account the microscopic selection 
rules. 



Here 
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Appendix A: Analytical results for the 
charged-phonon theory of IR response in ungapped 
bilayer graphene 

In this Appendix we provide the analyical derivation 
of the optical properties of the A mode in the absence of 
electrostatic bias between the layers, A = 0. To this aim 
it is convenient, for the simple tight-binding model with 
7o~7i considered, and in the linear expansion close to the 
K point, to employ the cylindric coordinates, where the 
Hamiltonian reads: 



/A/2 



H k = 



V o 



A/2 

7i 






7i 

-A/2 
et. 



-A/2/ 



(Al) 



where e k = tkC e = hv\k\, and 9k = arctan(fc y / k x ) . 

The Hamiltonian is diagonalized by the transforma- 
tion: 



where 



E,, 



H k = M k 1 H k M k = E k , 



[E lk \ 

E 2k 

E 3k 

V E 4k J 



(band labels according Fig. [Tp), 

M k = Ro^Fk, 



( 



> i8k 

10 

1 

e~ lB 




-Cfc Cfc s k 

Sk s k c k 

Sk ~s k c k 

-c k -c k s k 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



2sfcC fc 



(A9) 



In the diagonalized basis we have: 



M k " ! jk, 9 Mk 



/ -S k s g 
iC k C0 
CkSg 

-iSkCg 



-iC k Cg C k Sg 

-SkSg iSkCg 

-iSkCg SkSg 

-CkSg iCkCg 




where S k = 2s k c k , C k = c\ - s\, sg 
In a similar way we obtain: 



smy, eg — cost 



U 



k,A 




(All) 



/ -SkCg iCkSg CkCg 

-iCkSg -SkCg -iSkSg 

CkCg iSkSg Sk eg 

\ iSkSg -CkCg -iCkSg 




Note that the interband transitions between bands 1-4 
and 2-3 are missing in the electron-phonon matrix (jAllj) 
for the A mode, so that these interband transitions will 
not be operative in any response function involving such 
electron-phonon coupling, like for instance the mixed re- 
sponse function for the infrared activity as well as the 
phonon self-energy. 

The mixed response function reads in the diagonalized 
basis as 

Xju (iuim) = N S N V — ^2 Tr [jk,v9 iuj n + iw m ) 



'N ^ 

k . n 



(A13) 



where g(k, iuj n ) — l/[(ihuj n + fi)I — Hk], is the electronic 
Green's function in the diagonalized basis. 

Using now Eqs. ([A3]l , (|AT0|) . (|ATT|) . ([AT2]) , we can 
easily obtain: 



k,a,/3 



where 



n a/3 _ 



N s N v (j ki y) al3 (V k ^)/3 a , 



(A14) 



(A15) 
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and 



f{E k , a - fJt) - f {E k ,fi - n) 



Ei. 



E. 



fe,/9 



hibJr, 



(A16) 



Note that the function (iw m ) does not depend, in 
cylindrical coordinates, on the angle 9 k but only on 
the momentum modulus k — |k|. Writing ^ k — 

J 2nkdk J d9k/2ir, the quantity C°f k can be replaced 

thus with its average over the angle 9k, C"f k — > 



C 



(C Q/ 



a[3 



Using JXTTJ, ([A"12]) . it is now 



easy to see that C"§ k averages out, whereas \CJ A 

gevN s N v -fi / 'Ay/ '(hvk) 2 + 7^, which recovers Eq. (|41[) . 
The relative change of sign between different inter- 
band contributions stems from (ik,y)a/3(Mc,i/)/9o; = 



afi 



Appendix B: Charged-phonon theory for Raman 
spectroscopy 

In this Appendix we provide a brief derivation of the 
charged-phonon theory as applied to the case of Raman 
spectroscopy within the effective mass approximation. 
To this end it is useful to recall the tight-binding Hamil- 
tonian which can be written as 



H P = 



( A/2 7o /; 

7o/ P A/2 7l 

71 -A/2 7o/; 

V 7 o/p -A/2/ 



(Bl) 



where 



/ p = e-^ a/VI + 2e^ a/2%/5 cos(p y a/2). (B2) 



Close to the K=(47r/3a,0) point we can write p = 
K + k, and by expanding for small k we obtain 
([T)l, The Raman vertex 7 (0, </>') for a particular po- 
larization geometry can be now defined as 7 (0, <p') — 

a/N) E k , ff *Llk(& ^)*k,a where 



7k = (ei • V k )(e • V k )ff k 



= cos (j) cos (fj ,A f£ x 



+ cos 6 sin 



if ~ yx 



sin (j) sin 0' 7k a 
cos <^> sin (j) 1 ^ , 



(B3) 



where e; = (cos 0, sin 0) and e Q = (cos (j)' , sin </>') are the 
directions of the incoming and outcoming photon respec- 
tively, and 7 k J = d 2 Hk/dkidkj. Using Eq. (|Blj) we ob- 
tain, close to the K point, 



7g x (K) = - 7 f(K) 
7f(K) = 7f(K) = 
while, at the K' point, 

7 r(K') = -7f (K') 
7k y (K') = 7f(K') = 



-i{a v )/la\ 

= I(& x )/Aa 2 , 
i(a y )/4a 2 . 



(B4) 
(B5) 

(B6) 
(B7) 



We can write the Raman (RM) response function in 
the Matsubara imaginary time for generic polarization 
as: 

XRM(T,<t>,4>') = -(TrrtTtWMM'))- (B8) 

For unpolarizcd Raman scattering we have 
to average over <f> and <fi', Xrm(*w to ) = 
J d<f>/(2n)d^/(2Tv)xBM(Wm, $,<!>')■ We get then 

XRM(*w m ) = X7 xx 7 xx (i^m) + X^yy-iyy (iw m ) 



i—x,y 



where we used the shorthand notation x — y, y = x. 

For the electronic Raman scattering, taking also in 
account the different K and K' points, these four con- 
tributes are degenerate and we get on the real frequency 
axis: 



XrmM 



4Xn 



(BIO) 



Note that, since 7^ cx I{a x ) cx jk,x, this results explic- 
itly shows also that the electronic Raman background is 
directly related to the electronic optical conductivity. 

Let us focus now on the onset of phononic peaks within 
the framework of the Fano-Rice theory. For a generic 
case, for finite charge concentration n and finite gap A, 
we can write: 

i=x,y vv' 

+ E x^MA^Mx^^XBii) 

i— x,y 

where the label v specifies at the same time the A vs 
S phonon branch and the x vs y direction of the lattice 
displacement. The notation i denotes in addition i = y 
if i = x and i = x if i = y. 

Expanding explicitly the sums over i, i, v, and v' in 
Eq. (|B11[) . as well as the summation over the differ- 
ent K and K' points of the Brillouin zone, and taking 
into account the degeneracies j£ x = — 7 k v , 7 k y = 7 k x , 
-Ds^Sx = D SySy , D AwAi = D Ay A y , D Ax s x = E Ay s y , 
-C>SxA x = E)g yAy , after few straightforward steps, we end 
up with 

Axrm(w) = x 7 s(w) J D ss (w)xst 7 (w) 

+X 1 A(oj)D AA (uj)xAi~ f (u) 
+X~,s(u)D SA (uj)x A t~ / (uj) 
+X 7 AM£AsMxst 7 M, (B12) 

where, due to the degeneracies of the systems, the Ra- 
man vertex operators and the electron-phonon vertex op- 
erators have been chosen for simplicity in Eq. (|B12|) as 
7 = 7*» S = S x , A = A x . 
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